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ABSTRACT 

(— | 

Context. Due to its great mass and its rapid formation, Jupiter has played a crucial role in shaping the Solar System. 
The knowledge of its internal structure would strongly constrain the solar system formation mechanism. Seismology is 
Q ■ the most efficient way to probe directly the internal structure of giant planets. 

Aims. SYMPA is the first instrument dedicated to the observations of free oscillations of Jupiter. Principles and theo- 
^5 ' retical performance have been presented in paper I. This second paper describes the data processing method, the real 

instrumental performance and presents the first results of a Jovian observation run, lead in 2005 at Teide Observatory. 
Methods. SYMPA is a Fourier transform spectrometer which works at fixed optical path difference. It produces Doppler 
shift maps of the observed object. Velocity amplitude of Jupiter's oscillations is expected below 60 cm s _1 . 
Results. Despite light technical defects, the instrument demonstrated to work correctly, being limited only by photon 
noise, after a careful analysis. A noise level of about 12 cm s _1 has been reached on a 10-night observation run, with 
21% duty cycle, which is 5 time better than previous similar observations. However, no signal from Jupiter is clearly 
highlighted. 

Key words, planets and satellites: formation, Jupiter: oscillations, methods: observational, instrumentation: interferom- 
eter,techniques: spectroscopic 
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f^) . 1. Introduction their infrared detectors, not enough sensitive, did not detect 

^" ■ any signal. On the other hand, oscillations were sought in 

t j> ' Due to its great mass and its rapid formation, Jupiter has ve l oc ity measurements, obtained by Doppler spectrometry. 

k> I played a crucial role in shaping the Solar System. Two see- Schmider et al. (1991, hereafter S91) used sodium cell spec- 

' narios are generally proposed for the formation of giant trometer and Mosser et al. (1993, 2000, hereafter M93 and 

C3 ; planets: the nucleated instability (Safronov & Ruskol 1982) M00 ) the Fourier transform spectrometer FTS (at CFHT, 

and the gravitational instability models (Cameron 1978 and Hawaii) at fixed optical path difference. An excess of power 

Mayer et al. 2002). An efficient constraint on the formation has heen brought out in the spectrum at frequency range 

scenario would be given by measuring the total amount of [ 8 _ 2 ] mHZ) as well as the large separation of oscillation p- 

heavy elements inside Jupiter and the size of the planetary modes ar0lmd 140 ^Hz. Nevertheless, the oscillation modes 

core. Moreover, the knowledge of Jupiter's internal struc- have never been individually identified, hindering any con- 

ture would constrain the high pressure hydrogen equation of straint about the internal structure 

state, which is still inaccurate, and would particularly solve ' S YMPA is an instrument dedicated to Jovian oscilla- 
the question of the nature of the metallic-molecular phase tionS) which concept and performance have been described 
transition (e.g. Guillot et al. 2004). Gudkova & Zharkov in paper j ( Sc hmider et al. 2007). For the first time, a spe- 
(1999) showed that the observation of oscillation modes up cific i nstrumen t dedicated to Jovian oscillations was devel- 
to degree I = 25 would strongly constrain Jupiter's inter- Qped; inclu di ng imaging capability. Indeed, full disc obser- 
nal structure by exploring both the hydrogen plasma-phase vations do not permit to distinguish modes of degree higher 
transition and the supposed core level. than 3 More o V er, the broadening of the solar lines, due 
Attempts to observe Jovian oscillations have been t o the fast rotation of Jupiter, reduces drastically the sen- 
brought since the mid 1980's, thanks to different tech- gitivity of such measurements. The instrument, a Fourier 
niques. On one hand, Deming et al. (1989) have looked tachometer, is composed of a Mach-Zehnder interferome- 
for oscillation signature in thermal infrared. Unfortunately, t er which produces four images of the planet, in the visi- 
ble range corresponding to three Mg solar absorption lines 
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Fig. 1. Total intensity of each image. Nights 2 and 8 and 
the end of night 4 present a strongly and rapidly variable 
amount of photons, because of cloudy conditions. 
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Fig. 2. Spectral window of the observation campaign. The 
amplitude has been normalized to the spike's power. Only 
23% of the total energy remains in the central peak. 



quadrature, allows us to retrieve the phase of the incident 
light, which is related to the Doppler shift generated by the 
oscillations. 

Two instruments were built at Laboratoire Fizeau (Nice 
University). Three campaigns were lead simultaneously on 
two sites: in 2003, at San Pedro Martir (Mexico) and 
Calern (France) observatories; in 2004 and 2005, at San 
Pedro Martir and Teide (Canaries) observatories. 2003 cam- 
paign was mainly dedicated to technical commissioning. In 
Canaries, bad weather conditions have strongly limited the 
efficiency of the 2003 and 2004 campaigns, whereas the 2005 
campaign benefited of better conditions. 

In this paper we present the data processing, the real 
performance and the first results, obtained during the 2005 
run at Teide observatory. The processing of the San Pedro 
Martir data and the combined analysis of both observing 
2005 network campaign will be considered in a future work. 
After a short presentation of the observing conditions (Sect. 
[5]), Sect. [3] presents the process of the data analysis. We ex- 
pose in Sect. 2] all the steps for an accurate calibration of 
the data. Sect. [5] is devoted to the further of the data reduc- 
tion, for the measurement of velocity maps. The analysis of 
the time series get during the 2005 run at Tenerife is ex- 
posed in Sect. [6l Section [8] is devoted to conclusions and 
prospect ives. 



2. Observations 

Observations were conducted at Teide observatory 
(Canaries islands), with the 1.52-m Carlos Sanchez tele- 
scope, between march 31st and April 10th 2005. As it has 
been detailed in paper I, four images of Jupiter come out 
from the instrument. Optical parameters inside SYMPA's 
box are arranged such as the four 1.3-arcmin fields, cover 
128 pixels on the receptor (a DTA CCD, 1024 x 256 pixels). 
During the run, Jupiter at opposition presented a diame- 
ter of 48 arcsec, corresponding to 69 pixels on the CCD 
camera. 

Seven nights over ten days benefited of good weather 
conditions, yielding to a 21% duty cycle. Data quality was 
almost constant from one night to another, excepted for 
night 2 and 8, where clouds have reduced the incident 
flux. Observation conditions are summarized in Fig. [T] and 
Table [TJ The window function associated to the whole cam- 
paign was very strong, since we consider only Canaries data. 
Therefore, in the power spectrum, the amplitude of a single 
spike has been divided by a factor 4, since its power gets 



diluted in high side lobes (Fig. [5]). The total flux was ex- 
pected to be about 2.4 10 9 photons per 6-s exposure (paper 
I). However, its mean value along the run is about 2.2 10 9 
photons per exposure (Fig. [1} . This discrepancy with the 
estimated flux introduces a factor 0.95 in the ratio to noise 
value. In paper I, a noise level of about 4 cm s _1 was ex- 
pected for a 16- nights observation campaign, with 50% duty 
cycle. Instead of such a performance, by considering only 
Canaries data, the noise level is therefore expected at 10 
cm s _1 . 

3. Data processing strategy 

3.1. Four interferograms in phase quadrature 

SYMPA's instrumental principles are fully explained in pa- 
per I and are summarized in Fig. ([3]). The four output 
beams can be described in the detector coordinates {x, y) 
by the following approximations : 
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where Iq is the continuum component of the incident light, 
that is to say the Jovian figure, 7 the fringe contrast and 
4>{x, y) the incident wave phase map : 

<f>(x,y) = 2™ Q b(x,y)(\ + ^) (5) 

where o~q is the central wavenumber of the input filter, 
A(x,y) is the optical path difference (OPD) and «d and c 
are the Doppler and the light velocities. The Doppler shift 
of the solar Mg lines comes from the combination of the rel- 
ative motion of Jupiter to the Sun Wj/s, relative motion of 
the observer to Jupiter «e/j + VE, ro t (distance between the 
two planets and Earth's rotation), Jupiter's rotation i>j iro t 
and, finally, the oscillations v osc . In the following, we write 
the Doppler velocity as the sum of : 



2 (vj/S +VE/J + «E,rot + «J,rot + Vosc) 



(0) 



The factor 2 is due to the fact that the Doppler effect gets 
doubled after reflection on Jupiter's atmosphere. The orders 
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Table 1. 2005 run at Teide Observatory 



Starting dcit6 


Ending deiL.6 


Duration 


JVIean sampling 


Number of 


Nlean Intensity 


Standard deviation 


m d h m s (UT) 


m d h m s (UT) 


h m s 


s 


acquisitions 


photons image -1 


photons image -1 


Apr-02, 23:48:14 


Apr-03, 07:03:40 


7:09:35 


5.45 


4727 


1.89 10 a 


0.67 10 s 


Apr-03, 23:30:14 


Apr-04, 07:02:29 


7:31:59 


7.13 


3792 


1.79 10 9 


5.27 10 8 


Apr-04, 23:16:12 


Apr-05, 06:58:17 


7:39:30 


6.75 


4080 


2.50 10 9 


0.76 10 8 


Apr-05, 23:25:29 


Apr-06, 06:50:23 


7:22:12 


6.68 


3971 


2.44 10 9 


3.80 10 8 


Apr-09, 23:20:17 


Apr-10, 06:31:59 


6:37:48 


6.28 


3799 


2.06 10 9 


3.62 10 8 


Apr-10, 22:51:11 


Apr-11, 06:33:11 


7:41:48 


6.33 


4376 


2.38 10 9 


1.14 10 8 


Apr-11, 21:09:35 


Apr-12, 06:30:17 


7:20:42 


6.83 


3872 


2.36 10 9 


0.97 10 8 




Mach-Zehnder interferometer 



Fig. 3. Schematic view of the SYMPA instrument. The 
incident light coming from the 1.5-m telescope, passes 
through a 120-mm collimator and the 5 nm bandwidth in- 
terference filter. The optical path difference A occurs inside 
the Mach-Zehnder prism; it is a function of the heights H 
and h, refraction index N and n and incidence angles into 
the prism / and i : A = 2(H N cos / — hn cos i. The wollas- 
ton polarizer separates each output from the interferometric 
device into two separated beams. In total, the instrument 
produces on the camera four images of the same field, sep- 
arated by 7r/2 in phase (see paper I). 



of magnitude of these different terms are presented in table 

m 

From Eqt. the phase map appears to be the sum of 
two contributions : a "motionless" term, 2naoA(x,y), and 
a velocity term, 47rcroA(a;, y)vr>/c. The art of extracting the 
oscillation signal resides in the ability of eliminating step 
by step the motionless fringes and the "spurious" veloc- 
ity fields wj/s, «e/j, u E,rot and wj, ro t- All the steps of the 
schematic process presented below are detailed in Sect. 4 
and 5. 



Table 2. Orders of magnitude of different Doppler shifts 
during the 2005 run at Teide Observatory. 





Velocity (m s 1 ) 


Jupiter- Sun 
Jupiter-Earth 
Earth rotation 
Jupiter rotation 
Oscillations 


0.4 

3142 
[-409, 409] 
[-12570,12570] 
< 0.6 



3.2. Extracting the oscillation signal 

Let us consider a quadruplet of interfering images on the 
detector's field. The differences between the two couples of 
images, which are in phase opposition, allow to cancel the 
continuous component of the interferograms Iq, in order 
to keep the interfering patterns. We write U and V the 
normalized interfering patterns : 

U = — ^ ^ oc 7 cos 6 (7) 

h+h 

V = — cx 7sin0 (8) 

h+h 

where we let go of (x, y) dependence in order to simplify 
notations. The incident wave phase is retrieved by taking 
the argument of the complex interferogram: 

Z = U + iV oc 7e'* (9) 

The data processing expands in three main steps: correction 
of the motionless fringes, elimination of Jupiter rotation 
and elimination of the relative motion of the observer to the 
target and of the target to the Sun. Mathematically, it lies 
in creating successively four complex interferograms, asso- 
ciated to each step of the data processing: Zq for motionless 
fringes, Zj^ ot and i?E,rot for Jupiter's and Earth's rotation, 
Z-e/j for the relative motion of the Earth to Jupiter and 
Zj/s for Jupiter's motion with respect to the Sun. Then, 
the rough Jupiter complex interferogram Zj up comes decon- 
volved from additive signals: 

Zjfi a t = Zjup x Zq x Zj Iot Zft mt x Z^/j x Zj/ S (10) 
cx exp (47ra A^) (11) 

where the asterisk indicates the complex conjugation. The 
resulting interferogram is so-called Zj^at because of the 
appearance of its argument, i.e. the velocity map. Indeed, 
since oscillation amplitude is expected to be lower than 0.6 
m s -1 , it is absolutely impossible to see directly oscilla- 
tion modes in a single phase map, where mean noise level 
is expected to be about 900 m s -1 per pixel (Paper I). 
Oscillations might be picked out only in the spectrum of 
long time series. 

Note that Doppler shifts due both to Earth relative 
motion to Jupiter and Jupiter relative motion to the Sun 
are uniform across the field. On the contrary, the Jovian 
Doppler component is very sensitive to Jovian rotation 
since its value varie from —12.57 km s -1 to 12.57 km s -1 
from east to west on Jovian equator. In terms of velocity, it 
implies that the error in positioning the phaser Zj.rot must 
be smaller than the photon noise. Therefore, it supposes 



4 



Gaulme et al.: SYMPA: II. Real performance and first results 





20 40 60 80 100 120 



Fig. 4. Simulation of intcrferograms along the data process- 
ing chain. Top: motionless interferograms (Uo, Vq). Middle: 
Jovian interferogram (C/jup, Vjup) when Jupiter is inclined of 
—72° with respect to vertical. The interference pattern is 
mainly due to the coupling between motionless fringes and 
Jovian rotation (Eqt.[6|). Bottom: Jovian interferograms de- 
convoluted from motionless contribution, (U TO t, V T0 t), i.e. 
fringe pattern associated to Jovian rotation. Note that 
fringes associated to solid rotation present velocity iso- 
values parallel to the rotation axis. Differential wind profile 
has not been introduced. 



to know Jupiter's position better than 1/20 th of pixel. 
Figure [4] presents simulations of the motionless interfero- 
grams, Jovian interferograms and Jovian rotation interfer- 
ing pattern. 

4. Data calibration 

4.1. Pre-processing operations 

Pre-processing consists in cleaning each quadruplet of 
Jovian image (7i, I2, I3, 14) in order to create couples of 
interferograms (U, V). This implies three main operations. 
First, the camera dark current contribution has to be 
subtracted, by using offset images. Second, the inhomo- 
geneities of the single pixel responses to light intensity (pho- 
ton/electron gain) have to be compensated by dividing each 
image by a fiat field image. Third, the construction of the 
Jupiter phase map, using the argument of the complex im- 
age created by the difference of two couples of images, re- 
quires that all the four images overlap one to each other. 




Fig. 5. Overlapping of grids 1 and 3, before then after dis- 
tortion rectification. The position of the horizontal and ver- 
tical lines are determined by fitting the zeros of the grid 
image derivatives, along both directions x and y. Second 
order polynomials arc sufficient, since rotation and barrel 
distortion do not require higher orders. Then, each intersec- 
tion point coordinates are calculated by solving numerically 
the four order equation coming from the combination of the 
vertical and the horizontal line fits. The mean distance be- 
tween intersections points is equal to 5.32 % pixel -1 



If the two first points are easy to realize, the last one is 
pretty delicate because of the required accuracy of about 
1/20 pixel. 

As for every optical system including lenses and prisms, 
field distortion is unavoidable. Although Jupiter was po- 
sitioned as close as possible to the optical axis, its large 
diameter involves differential distortions between the four 
images, making the overlapping impossible. The whole dis- 
tortion effect is supposed to be composed of only transla- 
tion, rotation and barrel distortion. This problem has been 
anticipated by putting a regular grid, engraved on a glass 
slice, at the instrument focus. The grid intersection posi- 
tions are used to characterize the distortion (Fig. [5]) . 

Let us consider one image among the quadruplet 
Ii(x,y), i £ [1,4], where (x,y) are the detector coordinates 
(CCD pixels). Because of optical distortion, the image value 
on the (xk,yk) point, k 6 [1,128], actually corresponds to 
the (x' k ,y' k ) point. The repositioning algorithm must pro- 
duce for each image I{ an image I[, defined in the detector's 
coordinates, as: 



I[{x,y) = Ii(x',y') 



(12) 



The distorted coordinates (x 1 , y') are related to the regular 
detector coordinates by the relations: 



f(x,y) 



and 



(13) 



where / and g are polynomials expressed as J2k Ck&y k ~^ 
where k is the polynomial order and j < k. The polyno- 
mial coefficients Ck are obtained by minimizing the differ- 
ence between the intersection point coordinates of the grid 
associated to the considered image Ii and the intersection 
point coordinates of the regular grid. The knowledge of the 
set of distorted coordinates (x 1 , y') is used to build the new 
rectified image I'^x^y), by interpolation of Ii(x',y') upon 
the regular detector coordinates. 

Interpolation is realized with cubic method. In Fig. [5] we 
present a couple of 2-grid images before and after reposi- 
tioning. Note that the translation repositioning represents 
the zero order of the transformation. A way to evaluate 
the accuracy of the repositioning process is to apply it to 
the previously repositioned grids; they should overlap (Fig. 
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right). The mean distance between the repositioned in- 
tersection points from one grid to another belongs to the 
range [1/20,1/15] pixel, which almost fulfills the required 
accuracy. 

In order to respect the conservation of the flux, correct- 
ing operations have to be processed in the following order: 
first, subtraction of the offset to the image and to to the 
flat-field, then rectification of the offset-corrected image and 
flat-field. 



Iprc— processed 



[I(x,y,T)-0(x,y,T)} 



rectified 



[F(x,y,T) - 0(x,y,T)] t 



(14) 



cctificd 



where / indicates the considered image (e.g. Jupiter), O the 
offset image and F the flat-field image. T stands for the 
temperature of the camera. Note that assessing the four 
output intensities are not strictly equal, no photometric 
balance has to be performed since it is implicitly done by 
dividing each image by the flat field. 

4.2. Motionless fringes calibration 

The next step of data processing is the deconvolution of 
Jovian complex interferograms Zj up from motionless inter- 
ferogram Zq. Therefore, we have to characterize the mo- 
tionless phase term 2waoA(x,y) (Eqt.O. Furthermore, as 
it has been reported in paper I, the four output beams are 
not in perfect phase quadrature: the discrepancies between 
actual measurements and theoretical expectation are about 
e = 28° . This shift compared to quadrature has to be quan- 
tified precisely across the whole field (x,y). 

The knowledge of the optical parameters of both tele- 
scope and instrument permits to describe the motionless 
fringes (cf simulations on Fig. 0]), but not their imperfec- 
tions. Strong constraints on the motionless fringes come 
from solar light scattered by the telescope dome. Indeed, 
excluding the spectral Doppler shift, it presents the same 
spectrum as Jupiter, since 517-nm magnesium lines are so- 
lar reflected lines. Moreover such a process enables to en- 
lighten all the detector surface. Actually, the only difference 
between motionless fringes and "sky" fringes is the uni- 
form velocity field, introduced mainly by Earth rotation, 
much less by Sun-Earth distance variation and by unde- 
sired thermal effects (see Sect. H]). These parameters are 
taken into account thanks to the IMCCE ephemeride data 
base (www.imcce.fr). 

6-hour long "sky" shots were taken on April 2nd. In 
order to increase the signal to noise ratio of sky inter- 
ferograms, images are averaged over one minute intervals, 
which corresponds to three images. Figure ^ presents one 
of the 330 couples of 1-minute exposure interferograms. 
Fringes are well defined and fringe contrast (~ 0.4%) is 
smaller than expected (0.8 %). Unfortunately, it appears 
that the map of both interferograms arc not plane, but 
show a parabolic curvature. The origin of such a residual 
signal is explained by the fact that sky shots are taken dur- 
ing the day: some direct light enters into the instrument, 
and does not follow the correct optical path. 

A method has been developed to get round this unfit- 
ness. Let us take two couples of interferograms (U, V) and 
(U' , V'), taken at two different dates. The phase difference 
S<j) between them is due to terrestrial motion and temper- 
ature variation. Their expressions are given by : 

U = U + 7u cos(0) U' = U + 7u cos(0 + 5<jy) (15) 




50 100 

x (pixel) 



50 100 

x (pixel) 



Fig. 6. Top: interference pattern (U, V) obtained with solar 
light scattered by the telescope dome, over one minute (sum 
of 3 images). Bottom: cut along x-axis of U and V, at y = 
55 pixel. Note that both interferograms present a varying 
contrast between 0.3 and 0.5 % and a significant unflatness. 
Moreover, they are not centered around 0, that makes the 
phase <fi = arg(C7 + iV) impossible to recover. 
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Fig. 7. Cut of interference patterns (At/, AV) along x-axis. 
After correction of the low spatial frequency term, fringe 
contrast still varies strongly, from 0.2 % to 0.6%. 



V = V + 7v sin(0 + e) V = V a + 7v sin(</» + 5<f> + e)(16) 

where (Uq,Vo) describe the unflatness, and (7u,7v) the 
fringe contrasts. The subtraction of these signals eliminates 
the unflatness terms and brings to the two interfering pat- 
terns : 



AU = V -U = -7u 



AV 



v'-v 



7v 



2 sm — 
2 



2 sin 



sin(y>) 



cos(ip + e) 



(17) 



(18) 



where ip is the mean phase ip = <f> + 5cj)/2. Figure [7] shows 
cuts of (AU, AV) along the x-axis; the unflatness problem 
has been corrected. The fringe contrast presents optimized 
variations across both fields, from 0.2% to 0.6%. 

The motionless complex interferograms Zq = Uq + iVq 
is obtained by fitting both quadrature shift and phase on 
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Fig. 8. Left: plot of the interference pattern AU as a func- 
tion of the other AV. The spiral structure comes from the 
fact that both fringe contrast and ellipticity vary across the 
field. Right: plot of the interference pattern AU' as a func- 
tion of the other AV' . Points are distributed along a circle 
(solid line, white), which indicate that the two interfero- 
grams are in phase quadrature. In terms of phase accuracy, 
the standard deviation of these points is 5.3°. 



the couple (At/, AV). First, the plot of one of the inter- 
ference pattern as a function of the other highlights the 
phase quadrature imperfections (Fig. [8]). In case of perfect 
quadrature and uniform contrast in both interferograms, 
points would be distributed around a circle centered on 0. 
As a result, points are distributed along an ellipse, whose 
radius varies strongly on detector field. The phase shift e 
with respect to quadrature is related to the ellipse param- 
eters (A, B) by: 



arcsm 



A 2 -B 2 
A 2 + B 2 



(19) 



Since ellipticity and amplitude vary across the field (Fig. [51 
left), a further correction has to be applied. Therefore, the 
field (x, y) is divided in 10-pixel large squares, in which all 
these parameters are considered as uniform. Ellipse axis are 
estimated by least square fitting. The resulting parameters 
estimate (^4, B) 7 obtained for each sub-region, are interpo- 
lated for each pixel, by fitting their values by a 4th order 
polynomial. Thereafter, thanks to A, B and e, the ampli- 
tude of both interference patterns (AU, AV) is normalized 
to 1 and the phase shift is set to 90°, by the operation: 



, _ At/norm COs(g/2) - AKorm sin(g/2) 

cos(e/2) 

, AKorm cos(e/2) - AC/ norm sin(e/2) 



AV 



cos(e/2) 



(20) 
(21) 



where the subscript "norm" indicates that amplitudes have 
been normalized to 1. These new variables are now in the 
required quadrature. Figure [5] shows the plot of the in- 
terference pattern U as a function of the other V; phase 
quadrature is reached. At last, the motionless phase </>instru 
is obtained by fitting with a 4th order polynomial the ar- 
gument of the complex interferogram Z s y y = AU' + iAV' . 
The resulting phase standard deviation is about 5.3°. 

The motionless interferogram used in the following to 
process Jovian data is simply built as: 



Z 



o^instru 



(22) 




Indeed, since only the phase term matters, the fringe con- 
trast is set equal to 1 in the whole field. 



200 250 300 

x (pixel) 

Fig. 9. Top: quadruplet of Jupiter images (Ji, 1%, I3, 14), 
taken during night of April 2nd 2005 at Teide observatory. 
Bottom: cut of the top images. The instrument exhibits a 
global polarization effect, resulting in a better transmission 
of one channel with respect to the other. 



5. From four Jupiter images to a velocity map 

Processing the data lies in turning each quadruplet of 
Jovian images into a calibrated radial velocity map. The 
first step, developed in the previous section, gives clean 
fringes on both interference patterns (Uj up , Vj up ), in order 
to create the complex Jovian interferogram Zj up . Second, 
motionless phase is deconvoluted to Zj up with the help of 
Zq (Eat. [2"2")l. Third, Jovian rotation and other uniform ve- 
locity drifts have to be deconvoluted, in order to extract 
the velocity map. 



5.1. Cleaning Jovian fringes 

Let us consider a quadruplet of 4 pre-processed Jovian im- 
ages. Since they have been corrected from flat field and op- 
tical distortions, each couple of images {Ii,h) and (l2,h) 
present the same intensity level and shall overlap (Fig. [5]). 
Then, the two interfering patterns {Uj up , Vj up ) are created 
following Eqts. and [5] Moreover, they are set in phase 
quadrature following relation (|20[) and (f2"Tj) . As it can be 
seen in Figs. (TIT))) interference fringes appear, but, as for the 
"sky" interferograms, the fringes do not oscillate around 0, 
but around a distorted surface. These features are photo- 
metric residues, which have not disappeared with opera- 
tions (0 and ©. However, fringes have to be repositioned 
around a flat surface, centered around 0, in order to follow 
the phase of the Doppler signal. 

The best way to separate the photometric noise and 
spurious signal resides in filtering in the spatial frequency 
domain. The two dimension fast Fourier transform (FFT) is 
applied to the complex interferogram Zj up . In order to avoid 
spectral leakage due to the finite size of the image, we ap- 
ply an oversampling on the data Z ]np , by a factor 2, before 
passing to the Fourier domain. Besides, Jupiter's bound- 
ary is apodized with a cos 3 function, to avoid Airy-like re- 
bounds due to Fourier transform. In Fig. 111! we present the 
Fourier transform modulus. The central region contains the 
low frequency information (mean value, slow distortions). 
The horizontal line corresponds to the interference pattern; 
the horizontal structure comes from the fact that optical 
path difference varies essentially with ^-coordinates. The 
inclined alignment of stains constitutes the photometric ef- 
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Fig. 10. Rough interferograms (Uj up , Vj up ) obtained follow- 
ing Eqts.[7]and[Hl The interference pattern is clearly visible, 
but some photometric residuals of same order of magnitude 
as fringe's contrast, related to Jovian bands and zones, still 



feet. Indeed, it corresponds to the spatial frequencies along 
the rotation axis of Jupiter, which is inclined of about —72° 
in the detector's field. These features are actually remnants 
of cloud zones and belts. 

Wiener filtering stipulates that signal can be filtered 
out from noise in the Fourier domain if they were clearly 
differentiable. As it can be seen in Fig. [TTJ fringes occupy 
distinct positions from photometric residues. Actually, for 
high degree modes, a coupling between photometric noise 
and oscillation signal still exists, since high degree spherical 
harmonics extend largely in the Fourier domain. Therefore, 
part of the energy in the largest spatial frequencies may 
be filtered during this operation. Nevertheless, higher fre- 
quency mode information will be still available, since most 
of its signal is not cancelled by filtering operations, but may 
suffer from amplitude estimate uncertainty. 

In order to limit the damage generated by filtering the 
noise, the spatial filter is as smooth as possible. It consists 
of an ellipse, which major axis inclined of —72° with respect 
to vertical, which includes only the photometric stains and 
the central region. Besides, as for Jupiter, filter's bound- 
ary is apodized to avoid rebounds by applying the inverse 
Fourier transform, when getting back to the image plane. 
The resulting couple of interference patterns, after spatial 
filtering is presented in Fig (fPCT)) . Now, fringes are centered 
around the value. Hence, requirements for deconvoluting 
the motionless phase are achieved. 

5.2. Jupiter velocity map 

The motionless fringe deconvolution is realized by the fol- 
lowing operation: 



•2j,rot — -^jup x Zq 



(23) 



The resulting interference pattern presents fringes parallel 
to the planetary rotation axis (Fig.[T3"|). Indeed, the projec- 
tion of Jupiter's velocity field v = flB,(L), where fl is the 
angular velocity and L the latitude, towards the observer 
reduces to v = fix, where x is the abscissa along Jupiter's 
equator. At zero order, Jovian rotation can be considered 
as solid rotation since differential rotation with respect to 
solid-body rotation is about 1% on the equator. Thus, ft 
is almost uniform on the Jovian disk, and the phase of the 
complex interferogram presents iso-values along the rota- 
tion axis. Moreover, note that since noise level is about 900 
ms _1 per pixel, differential rotation is definitely invisible 
on a single image. 



photometric residues 




fringes 



k x (pixel" 1 ) 

Fig. 11. Modulus of Fourier transform of the complex in- 
terferogram Zj up = [/j U p + iVj U p. Photometric residuals are 
located along the direction indicated by the two upper ar- 
rows. This direction fits with Jupiter's inclination in the 
detector's field. Main fringe information relies in the hori- 
zontal line, under the central region. 
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Fig. 12. The same couple of interferograms after filtering 
of the photometric residuals. Note that fringes oscillate 
around 0, and that contrast is about 0.6%, which is very 
close to expected values. 
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Fig. 13. Interference pattern after deconvolution of the mo- 
tionless component. Fringes present iso-values along the ro- 
tation axis of Jupiter. The white circle indicates Jupiter's 
size before resizing. 



Hence, a complex phaser which reproduces the solid ro- 
tation is applied to each Jovian complex interferogram; it is 
defined by Z so iidrot = exp(47rerAv rot /c), where v rot = 2tt/T 
and where T = 9 h 55 m 30 s (system III) is the mean rota- 
tion period. The main difficulty of such a process lays on the 
accuracy of the estimate of Jupiter's position on the detec- 
tor, because 1 pixel corresponds to 350 ms -1 . Two meth- 
ods have been envisioned to make the center of the solid 
rotation phaser ^soiidrot overlap on the center of Jupiter's 
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interferogram Zj, r ot- In both cases, a threshold is applied 
to each image, in order to get rid of spurious photometric 
signals, such as Jovian satellites of terrestrial atmospheric 
light scattering. The "barycentcr" method consists in tak- 
ing the coordinates of the barycenter of the total photo- 
metric image / = Yli=i^i- The "interspectrum" method 
consists in determining the relative distance between two 
images, taken at different dates, by measuring the phase <& 
of the interspectrum of the two images. Indeed, the phase of 
the interspectrum of the couple of images (Ji,^) is defined 
as: 



$ = aig{T{h) x T{h)*} 
oc O2 - xi) + [y 2 - yi) 



(24) 
(25) 



where (a?i,yi) and (2:2,2/2) are the coordinates of the cen- 
ter of Jupiter and J- and T* indicate the Fourier and in- 
verse Fourier transformations. The second method has been 
preferred because the barycenter estimate is too sensitive 
to high spatial frequency photometric details, which vary 
along the night, as cloud features or satellite transits. On 
the contrary, the interspectrum method is sensitive only to 
low spatial frequency. 

After deconvolution of Jovian mean rotation (i.e. solid 
body approximation), interferograms become flat since the 
remaining phase (f> = Ana (y 3 / S + v E /j + w E ,rot + v osc ) is 
uniform across the field (Fig. Q3]). The subtraction of Wj/s 
and fE/j is exposed in the next section, since they do not 
change the noise to signal ratio of the phase map. Therefore, 
the velocity map is retrieved thank to Eqt. 6 in paper I: 



v = u arg{Z flat } 



(26) 



with f — f km s _1 . Note that because of apodization cre- 
ated by spatial filtering of photometric residuals, the entire 
phase map is not exploitable. A part of the external region 
is cut down, whose proportion is a function of the apodiza- 
tion strength. Here, an external ring representing 1/4 of the 
Jovian radius is taken away (Fig. [TJ|. As a consequence, 
the flux is reduced by about of the partial elimination of 
Jupiter's is the loss of 43 % of the photons, which make the 
expected noise level go up around 10.5 cm s . However, 
the performance decrease is limited by the low weight of 
external regions in the Doppler signal. 

The standard deviation of velocity across the Jovian 
disk is about 890 ms^ 1 per pixel, that is to say 18.9 ms" 1 
when integrating the 2200 pixel of the resized Jovian disk. 
This performance matches with expectations. If photon 
noise is reached, a 7-hour night integration with 6-s sam- 
pling yields a noise level as small as 29 cm s . 

6. Temporal analysis over one night 

In the previous section, we have exposed the data process- 
ing method for a single Jupiter quadruplet, for the extrac- 
tion of the Doppler signal. In this section, we analyse the 
time series of the median phase extracted for each Jovian 
quadruplet. Then, we identify the different sources of noise 
and present correction methods. 



6.1. Phase global behavior 

The study of median phase along a single night allows us 
to evaluate the mean noise level and to identify spurious 
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Fig. 14. Left: Jupiter phase map, after deconvolution of 
motionless phase and Jovian rotation. Terrestrial motion 
and temperature variations have not been taken into ac- 
count, since they only shift the whole phase by a uniform 
value. Right: complex diagram of the same interferogram; 
Vfi a t is plotted as a function of £/fl a t- The white circle in- 
dicates the standard deviation of the phase of the interfer- 
ogram, which value is 8.910 -4 rad. In terms of velocity, it 
corresponds 890 ms -1 per pixel, that is to say 18.9 ms -1 
on the whole disk (2200 pixels are still available on Jupiter 
after resizing). 

signals. Median phases are extracted from Jupiter's flat in- 
terferograms Zflat as follow: 



V m 



mod 



(27) 



where the subscript "med" indicates the median value of 
the considered variable. This estimate of Jupiter median 
phase has been preferred to the direct median of the phase 
map, in order to avoid noise coming from 27r jumps (see 
Fig. [T4"|) . In Fig. [T51 median phase along night 9 and its 
power spectrum have been plotted. 

First, according to data processing chain (Eqts. [5l [TU] 
and Hip , phase measurement is proportional to the red- 
shift of spectral lines, it appears that the signal is over- 
whelmed by a strong low frequency noise. Indeed, along 
the night, the velocity variation should be dominated by 
the Earth rotation (409 m/s amplitude at Teide observa- 
tory), so the phase shall increase of about 0.41 rad instead 
of decreasing of about 1.12 rad. This implies that an unex- 
pected stronger low frequency drift dominates. Besides this 
low frequency noise, the power spectrum highlights a rapid 
oscillation noise source around 6-mHz frequency. 

Anyway, the standard deviation of median phase is 
around 27 m s _1 per image; images were taken every 6 
s. By supposing an only photon noise origin, the noise level 
reduces to 15 cm s _1 for 7 nights of 7 hours, which is 1.3 time 
worse than expected from previous section. In the following 
subsections, we interpret the noise origins and describe the 
method which has been used to get rid of their main effects. 



6.2. Temperature effect 

Paper I has expose that the interferometer is made of two 
pieces of two different glasses, specially chosen as to com- 
pensate the index variations and the dilatation, in order 
to have a stable OPD. In particular, a temperature varia- 
tion of about 1°C should have no effect on OPD in 14°C 
environment. 

The median of Jupiter's median phase over each night 
has been plotted as a function of the corresponding me- 
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Fig. 15. Top. Median phase along night 9: measured phase 
(full line) versus theoretical phase drift due to Earth rota- 
tion and relative motion of Jupiter toward the Sun and the 
Earth (dashed line). Bottom. Amplitude spectrum of the 
measured phase as a function of the time frequency. Two 
spike forests blow out the mean noise level at less than 0.5 
mHz and at 6 mHz. Note that mean noise becomes totally 
flat in the frequency range above 20 mHz, which means that 
photon noise is achieved, around 41 cm s . 



dian temperature fFig. I16[) . Phase and temperature appear 
clearly correlated. A variation of 1°C introduces a phase 
shift of 0.33 rad (i.e. 330 m s -1 ), which is much above 
than expected (between —60 m s _1 and 30 m s _1 for tem- 
perature between and 20° C; see paper I). In fact, such 
a thermal effect is retrieved when taking into account the 
error bars upon dilatation coefficient and refractive index 
thermal dependance (± 10 %). 

Unfortunately, the correlation coefficient reported in 
Fig. [16] does not allow us to correct Jovian phase map from 
its thermal dependence. Indeed, the comparison of phase 
and temperature over one night does not match correctly: a 
several hour delay is still present (Fig.[T7|). This discrepancy 
is due to the fact that temperature measurements do not 
correspond to the Mach-Zehnder prism, but to the metallic 
box which holds it on. Because of a greater thermal iner- 
tia, glass actual temperature is time-shifted with respect to 
metal temperature. Therefore, a low frequency filter is ap- 
plied to data in order to reduce the noise in the frequency 
range below 0.2 mHz. 

Beyond a mean OPD drift, temperature variations gen- 
erate differential OPD variations with respect to light in- 
cidence angle onto the prism. Since Jovian interferograms 
get deconvoluted from motionless fringes with the help of a 
constant interferogram fEqt. l2"3"|) . the differential OPD vari- 
ations introduce a slowly varying inclined surface in Jovian 
flat phase maps (Fig. [T8]) . This surface is bent along the 
x-axis because OPD is only x-dependent. Such an inclined 
surface yields a high frequency perturbation correlated to 
Jupiter's position in the observed field. This noise is re- 
duced by fitting the spurious surface on each phase map, 
by a plane inclined only with respect to x-axis. Thereafter, 
Jupiter's phase map are set flat by using a smoothed esti- 
mate of the fitted parameters. 



) = 0.33 T - 1 .77 
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Median measured temperature (°C) 

Fig. 16. Phase-temperature correlation. For each night, we 
have plotted the median of the median phase as a func- 
tion of the median temperature over the entire night. The 
straight line represents the least square estimate of the cor- 
relation between phase and temperature. 
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Fig. 17. Jovian median phase, corrected of Earth's rotation 
component, over night 10 (full noisy line) and temperature 
converted to phase with the previously fitted correlation 
coefficient (full smooth line) . Note that measured phase fol- 
lows the temperature with a several hour delay. 
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Fig. 18. Top: phase map averaged along y-axis as a func- 
tion of the x-axis. A negative slope is clearly visible. 
Bottom: root mean square estimate of the slope of the in- 
clined plane along night 10. 



6.3. Guiding noise 

The so-called "guiding noise" is the perturbation related 
to Jupiter's position, which should not occur if Jupiter 
were at fixed position. Theoretically, SYMPA's Doppler 
velocity measurements are not sensitive to Jupiter's po- 
sition onto the field, because after motionless fringe de- 
convolution (Eqt. [23]) nothing should depend on the co- 
ordinates. However, the comparison between spectra of 
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Fig. 19. Top. Jupiter's coordinates (x,y) in the CCD field 
over night 9. Note that y varies much more than x coordi- 
nate because Jupiter is rotated about —72° on the image 
and because most of the telescope guiding problems occur 
with the right ascension. Bottom. Amplitude spectrum of 
both coordinates. They present strong components below 6 
mHz and significative spikes up to 18 mHz. 



Jupiter's phase and coordinates (x, y) (Figs. [T5] and 
shows a strong correlation at 6 mHz. In fact, many other 
correlated signatures occur at many frequencies. 

This guiding noise has two main origins. First, the 
temperature high frequency noise cannot be totally can- 
celled because of the inclined surface fit inaccuracy. Second, 
Jupiter's position is estimated around 1/15 of pixel. 
Therefore, the deconvolution of Jovian rotation introduces 
a spurious signal correlated to Jovian position. Thus, the 
guiding noise is a combination of these two effects, for which 
it is hard to determine which of them predominates. 

The guiding noise is strongly reduced by removing 
the intercorrelation of Jupiter's phase and coordinates. 
The decorrelated phase with respect to both coordinated 
<^decorr/x/y is obtained as follow: 



</>dccorr/x = T - I x )F{<t>}} 

?Wcorr/x/y = ? M (1 — ly ^{^decorr/x}} 



(28) 
(29) 



where 0decorr/x indicates the phase decorrelated with re- 
spect to x, and I x and I y the interspectra relative to x and 
y, defined by: 



Iy = 



2 arcsin(jR{J r (^)J : -*(x)}) 

2 arcsin (^{^(^decorr/x)-? 7 * "(x)}) 

ecorr/x 



(30) 
(31) 



Note that we apply a 75% threshold upon I x and 7 y in order 
to decorrelate only significant guiding noise spikes. 

Jupiter's phase standard deviation drops down from 27 
m s _1 per image to 21 m s _1 , which is very close to the 
expected photon noise level (19 m s _1 per image). Hence, 
by supposing the remaining noise is only due to photons, 
the noise level reduces to 32-cm s _1 on a 7-hour integration 
and to 12 cm s^ 1 after 7 nights. The remaining not-photon 
noise is ignored. 



7. The search for oscillations: no evidence for a 
Jovian signal 

7.1. Decomposition into spherical harmonic base 

A stationary oscillation mode can be described as the sum 
of two spherical harmonics of degree I and order ±m. 
Jupiter velocity field related to p-modes follows such a de- 
scription. Therefore, the Doppler signature of the radial 
velocity field expands in a base made of the projected com- 
plex spherical harmonics towards the observer. 

The coefficients c™ associated to each spherical harmon- 
ics Y™ arc obtained as follow: 

m _ Epixcis^T 1 }^ ^ £ P ixoi g 3{>7"}x" ,„„, 

C l - mrvml + 1 rvrvml V* 1 ) 



Spixels ^07™} 



y . 

Z— 'PI 



pixels 



where v indicates the velocity map (Eat. [21)]). The normali- 
sation coefficient is done with respect to the actual number 
of pixel after resizing. The c™ coefficient modulus are ex- 
pressed in m s _1 . Oscillation search is performed in the 
spectra of all spherical harmonics up to the degree I = 25. 

7.2. Regular temporal grid 

Velocity maps are extracted from mean flat interferograms, 
averaged within 30-s intervals, but calculated every 15 s. 
Such a process has two reasons. On one hand, the mean 
noise level inside velocity maps drops down of a factor v5 
since 30 s contains 5 images, which limits strongly the 2-7T 
jumps which appear when applying Eqt. 1261 On the other 
hand, it allows to use the fast Fourier transform (FFT) to 
calculate the power spectra. The procedure is of great inter- 
est because the search for modes up to 25 means 625 spectra 
of 27000 point in the time series. The spacing of data ev- 
ery 15 s imposes a cut-off frequency at 16.3 mHz, which 
is well beyond the expected p-modes (less than 3.5 mHz 
from Mosser 1995). The averaging within 30-s intervals is 
done in order to avoid spectrum leakage when calculating 
the Fourier transform. 

7.3. Power spectrum of modes (£ = 0, m = 0) and 

(£ = l,m= 0) 

We present the power spectrum of the time series corre- 
sponding to spherical harmonics Yq and Y®, for all the data 
of Canaries observation campaign. The choice of these two 
modes among 625, permits to present the two main types 
of spectra. The first is sensitive to the remaining guiding 
noise, whereas the second is much less sensitive. Indeed, 
since guiding noise is mainly due to right ascension control 
defects, the mode {£ — l,m = 0) (hereafter (1,0)) is less 
sensitive to these problems since north and south Jovian 
hemisphere compensate (Fig. |2T)]). 

Both spectra exhibit a flat noise level in the frequency 
range [1, 8] mHz, excepted around 6 mHz, where a guiding 
signal subsists. Beyond 8 mHz, the averaging over 30 s cuts- 
off the signal. Guiding signature is reduced from 4 m s -1 
to 1.3 m s _1 after decorrelation processes. The mean noise 
level is about 12.6 cm s _1 for mode (0,0) and 11 cm s _1 for 
mode (1,0), which squares with the last-estimated photon 
noise level. Note that such a performance has never been 
reached on Jupiter and proves that the instrument and the 
data processing chain works efficiently. 
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Fig. 20. Power spectrum of time series related to modes (0,0), top, and (1,0), bottom, after concatenation of data along 
the whole run. The mean noise level are respectively of about 12.6 cm s _1 and 11 cm s , what matches with the last 
estimate of the photon noise. The 6-mHz guiding spike still is about 1.3 m s _1 for the (0,0) mode, whereas it is only 
about 0.5 m s _1 for the (1,0) mode. 



As regards the comparison to previous observations of 
S91 and M93 and MOO, no excess of power is present be- 
tween in the [1, 2] mHz frequency range. Moreover, no large 
spacing v§ is highlighted, whose value is estimated around 
150 /iHz, and which was detected around 136 and 143 /iHz, 
respectively, by S91 and MOO. Such a dissension with pre- 
vious observations will be analyzed in detail in a future 
work. 

As regards the excess of power in the frequency range 
[0.3, 0.6] mHz, no indication for a Jovian origin can be fur- 
nished at this step of the data analysis. It could be a re- 
maining low frequency noise, related to temperature and 
position. A global analysis over all the modes up to degree 
£ — 25 is required to determine its origin and to high- 
light global signature as Jovian rotation frequency or 
frequency. 

7.4. Global analysis: no evidence for a Jovian signal 

As for helioseismology, simultaneous temporal and spatial 
frequency analysis may reveal the presence of significant in- 
formation lost among noisy spectra. In Fig. [21] we present 
the (£, v) and (m, v) diagrams. In the (£, v) diagram, the low 
frequency excess of power, picked out in the (0, 0) and (1, 0) 
power spectra, is confirmed for almost all degrees £ in the 
frequency range [0.3, 1] mHz. It does not exhibit an orga- 
nized structure. Moreover, in comparison with past obser- 
vations, no signal is distinguishable in the [1, 2] mHz range. 
On the other hand, in the (m, v) diagram, the same excess 
of power appears to be strongly structured. The energy is 
distributed along 2 main lines, symmetric with respect to 
abscissa, beginning at 0.5 mHz and ending at 1.2 mHz for 
m = ±25. Moreover, a second couple of lines, almost par- 
allel to the first ones, is still visible between frequency 0.7 
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Fig. 22. Velocity maps averaged over 5 minutes. A fringe 
structure, inclined of about Jupiter's inclination on the de- 
tector is still present. The coupling of these feature with 
the spherical harmonic masks may explain the observed 
features in the (to, v) diagram. 



and 1.6 mHz. The mean slope of the principal lines is about 
28 ^Hz, which corresponds to Jovian rotation frequency. 

A quick analysis of the (to, v) diagram gives some indi- 
cations. A simple guiding noise origin is excluded, because 
such a spurious signal contaminates all the eigenmodes at 
the same frequency, as for the 6 mHz spike (see Fig. l2"Tj) . 
Furthermore, guiding noise becomes negligible beyond or- 
der to = 15, since its effect compensates when applying 
high order spherical harmonic filtering to velocity maps. 
At last, if the slope of about 28 /xHz could be casual, it 
could indicate a Jovian origin to the observed signal. 
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Fig. 21. Time-spatial frequency analysis. Top, the (I, v) diagram represents the mean of all the power spectra at a given 
degree £, as a function of all the explored degree. Bottom, the (m, v) diagram represents the mean of all power spectra 
at a given order m, as a function of all explored orders. Note the guiding noise is still present at 6 mHz, within a vertical 
line in both graphics. The maximum amplitude is around 1.2 m s . 



However, the Jovian signal hypothesis may vanish by 
supposing a coupling between two spurious signals. First, 
it can be noticed, after summation of velocity maps over 
10-min, that a fringe structure still remains, which peri- 
odicity fits almost with mo = 16 (Fig. |2"2"| . Second, it can 
reasonably supposed that photometric signal has not been 
totally removed after the Fourier filtering step (see Sect. 
15. ip . Thus, a signal modulated by Jupiter's rotation prob- 
ably underlies inside velocity maps. The coupling of these 
two spurious signals introduces a modulation of the signal 
by a cos(mo$) factor, where <f> = Q r t comes from the pho- 
tometric remnants; f2 r = 2Trv r indicates the Jovian rotation 
frequency. Consequently, when applying the spherical har- 
monic search algorithm (Eqt. [32]), a coupling appears be- 
tween the cos(mo < i ) ) modulation of velocity maps and the 
cos(m$) associated to Y™. With this assessment, velocity 
maps are modulated in the following way: 

v = I'd cos(£l r mot) cos(£l r mi) (33) 

= -y {cos[fi rot (mo + m)t] + cos[ft ro t(m - m)t]} (34) 

Therefore, a linear dependence v — m v r + C appears, 
where the constant term C = niQV r is about 450 /iHz, which 
matches with the observed origin of the two lines. 

8. Conclusion and prospects 

8.1. Instrumental performance 

The aim of SYMPA instrument was the detection and mea- 
surements of acoustic modes on the giant planets of the so- 
lar system, with a previously unequalled sensitivity around 
4 cm s . Such a performance was estimated for a 16-day 



observation campaign with 50% duty cycle (see paper I). 
By choosing to process only the Canaries data, since the 
instrument used in San Pedro Martir Observatory presents 
some defects which make the data more difficult to process, 
the duty cycle is only 21% over 10 nights. In this case, the 
noise level is reevaluated at 10 cm s" 1 . By taking into ac- 
count the lack of photons by a factor 2, underlined in Sect. 
[5J the 1— a sensitivity decreases to 12 cm s _1 . After decor- 
relation of time series with respect to Jupiter's position on 
the detector, outside the low frequency range (< 0.8 mHz) 
and the 6-mHz spike, the power spectra are flat and the 
mean noise level reaches 12 cm s _1 . Such a noise level is 5 
time better than previous observations which were limited 
around 60 cm s _1 (MOO). 



However, the data processing has highlighted some de- 
fects, as the too strong optical path difference dependence 
to temperature and some difficulties, as the separation be- 
tween photometric and spectrometric information. The lat- 
ter make the extraction of Doppler velocity hard to practise. 
This is mainly due to an insufficient accuracy in rectifying 
the distortion, particularly about the photometric effect of 
the distortion (variable PSF upon the detector). This point 
suggests some instrumental modifications, as the time mod- 
ulation of OPD in order to modulate the phase of the fringe 
pattern of about w. It would allow us to replace the spatial 
subtraction between interferograms by a time subtraction. 
Moreover, the spectral information would be emphasized 
more easily by narrowing the entrance filter, which would 
affect the global sensitivity. This would increase the fringe 
contrast. 
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8.2. Jupiter 

Jupiter's seismological observation of S91, M93 and MOO 
have all exhibited an excess of power in the frequency range 
[1, 2] mHz, which matches with the theoretical expectations 
(e. g. Mosser et al. 1996). Moreover, they did not provide 
any spatial resolution. If both observation methods (sodium 
cell and Fourier transform spectrometry at fixed OPD) have 
presented an excess of power in the same frequency range, 
the large spacing i/ estimate was very noisy and differed 
quite sensitively with respect to the theoretical value of 153 
^Hz (Gudkova et al. 1995). 

Our observations do not present features close to p- 
modcs signature: the absence of the large spacing vq in 
power spectra and (£, v) diagram is significative. However, 
it is worth to notice that because of the resizing of Jovian 
velocity maps, the projected spherical harmonic base be- 
comes a quite degenerated base. Therefore, information 
from a given mode goes diluted into other mode spectra, 
which may diminish strongly the p-mode signature and the 
identification of specific related structures. A way of leaving 
degeneracy has to be developed. 

Moreover, Bercovici and Schubert (1987) roughly esti- 
mated Jupiter's oscillation velocity amplitude between a 
few cm s _1 and 1 m s _1 . Therefore it is not abnormal that 
oscillations are not enlightened with a 12-cm s _1 noise level, 
with a 21% duty cycle. This reduces the observed amplitude 
of any oscillations by a factor (0.2) 05 . In these conditions, 
only a 25-cm s _1 signal could be detected at 1— a level. 

8.3. Prospects 

SYMPA has demonstrated to work properly, after taking 
into account its technical defects (temperature dependance, 
field distortions, difference of intensity between the two po- 
larized outputs, lower sensitivity of the CCD). It has per- 
mitted to reach the best remote sensing velocity measure- 
ments upon giant planets. Some improvements are to be 
performed on the existing instrument, as a better thermal 
insulation and a slightly modified optical design in order 
to reduce the geometrical distortions. An alternative and 
more efficient solution is to rebuilt a prism, equipped with 
an OPD time modulation. 

From a most general point of view, this seismometer 
is a tachometer, which furnishes velocity maps, instead 
of point to point measures, which is the case of cchellc- 
spectrometers. Thus, it can be used to other kind of ob- 
servations, such as wind velocity measurements. SYMPA 
has tested an extra seismological application in November 
2007, by participating to a ground-based observation cam- 
paign organised in sustain to ESA Venus-Express probe, in 
order to characterize the lower mesosphere wind velocity. 

What future for giant planet seismology is the natural 
question arising at the end of the first run of the first project 
entirely dedicated to this topic. The main obstacles to such 
a kind of measurements are the temporal coverage and the 
atmospheric turbulences. Temporal coverage introduces a 
windowing effect which makes the signal amplitude drop 
down and the deconvolution of noise with respect to signal 
very hard because of the spreading of information into large 
frequency ranges. Atmospheric turbulences limit the spatial 
resolution and, above all, make the planet move in the field, 
which generates guiding noise. A significant improvement 
would arise with Antarctic observations: Schmider et al. 



(2005) showed than a 80% duty cycle can be reached for 
more than 3 months. However, Jupiter south hemisphere 
oppositions will last till 2009, and afterwards will not occur 
before 2018. 

The ideal opportunity would come from space measure- 
ments aboard an interplanetary spacecraft cruising towards 
Jupiter. With an optical design similar to SYMPA, a 10-cm 
entrance pupil observing for 2 months at a mean distance 
of 0.2 AU to Jupiter would allow to reduce the noise level 
to few mm s _1 and to increase the spatial resolution till de- 
gree I = 100. Such an instrumental concept has been pro- 
posed to the European project of mission to Europa and the 
Jupiter system, Laplace, which participates to the ESA's 
Cosmic Vision programme (Blanc et al. 2006). The con- 
cept, so-called ECHOES, is derived from the ground based 
instrument SYMPA and the space helio-seismometer MDI 
onboard on SOHO space Solar observatory (Scherrer et al. 
1995). It would produce simultaneously intensity maps of 
reflected solar light, velocity maps and polarimctric maps of 
the whole visible surface with a spatial resolution of about 
500 km at the surface of Jupiter. Also, such an instrument 
would be an efficient tool to measure the wind velocity of 
the upper troposphere and the lower stratosphere. Such a 
spatial project needs a feasibility study. On one hand, from 
an instrumental point of view. Are the specific requirements 
of seismology compatible with the interplanetary probe pro- 
gramme? Is the envisioned method is the most adequate 
to space measurements, or COROT-type photometric mea- 
surements (Baglin et al. 1998) should be preferred as sug- 
gested by Mosser et al. (2004) and Gaulme & Mosser(2005) 
? On the other hand, theoretical improvements have to be 
performed about the modes excitation mechanism, in order 
to evaluate the amplitudes of the expected oscillations. 
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